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In the heavy quark effective field theory of QCD, we analyze the order l/mq contributions to 
heavy to light vector decays. Light cone sum rule method is applied with including the effects of 
l/niQ order corrections. We then extract \Vub\ from B — * plu decay up to order of l/mq corrections. 
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I. INTRODUCTION 



o 

o . 

. Much effort has been devoted to discuss the heavy to hght hadron semileptonic decays. In particular, B it{p)Iv 

^ ' decays attracted the most interest [1-9] because they can be used to determine the quark mixing matrix element 

O ■ IKifeli a parameter of significance in particle physics. The heavy quark symmetry and relevant effective theory greatly 

' simplify the study of hadrons each of which containing a single heavy quark and any number of light quarks, and 
provide relations between different processes. This symmetry is applied to study B{D)^g-^ ■n{p,K,K*)lv decays in 

^ Refs. [10-12], where the finite heavy quark mass (mg) corrections are not considered. Ref. [13] extends the study on 
B ttIv decay up to the next to leading order of the heavy quark expansion. For a more complete knowledge of the 

■ magnitude of the finite mass corrections to heavy to light meson decays, and to the determination of \Vub\, one should 

J> [ also study the l/mg order corrections to semileptonic B decays to light vector mesons. 

. In this short letter we will apply the heavy quark effective field theory (HQEFT) developed in Refs. [14-18] to 

' analyze the l/mq corrections to the B plv decay. And the light cone sum rule method will be adopted to 

] numerically estimate the nonperturbative functions, i.e., the heavy to light vector form factors with including l/mq 

^ order corrections. In section II the l/mq order corrections are formulated in HQEFT framework. Section III devotes 

ly-j to evaluate wave functions using light cone sum rule method in HQEFT. And section IV is the numerical results and 

• discussion. 

I ■ II. B pLu DECAY IN HQEFT 

Oh! 

(-H \ The transition matrix element responsible to the B — s- plv decay is generally parameterized by form factors as 



rriB + mn 



c5 ; +1 ^^y^ (e* ■ (p + qW + ^^y^ 6^"^^4(p + qhPy, (2.1) 

■ ■ niB + mp TUB + mp 

where q = pb — p is the momentum carried by the lepton pair. 

In the framework of HQEFT [14,15], the QCD quantum field Q for heavy quark is decomposed into particle field 
and antiparticle field Q~ , so that the quark and antiquark fields are treated on the same footing in a symmetric 
way. The effective quark and antiquark fields in HQEFT are defined as 



Q± ^ ^rMQV-XQ± ^ e'^rnQVXp^Q± (2.2) 



Rt - P^g± (2.3) 



with V being an arbitrary four- velocity satisfying = 1, and P± = (1 ± f)/2 being the projection operators. 
defined above are actually the large components of the heavy quark and antiquark fields respectively. are the 
small components of the heavy quark and antiquark fields respectively. The quantum field in QCD Lagrangian can 
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be written as Q = (5+ + Q~ = + Q~ + + R~ , which contains aU large and small components of particle and 
antiparticle. The decomposition of Q is presented in detail in Refs. [14,15]. 

After the small components of particle and antiparticle fields being integrated out, QCD Lagrangian turns into 



with 



^Q,v - ^'q,v '' + ^Q,v ^ + ^'q,v ■* + ^Q,t^ (2-4) 



^Q,t^ = QtiPvQt, (2.5) 



where 



2mQ V 2mQ 

ip± = ip- iini ■ D, = -0+ i ^v'd, (2.7) 

which is treated as HQEFT in the case that the longitudinal and transverse residual momenta, i.e. the operators 
iv ■ D and are at the same order of power counting in 1/ mg expansions. 

The heavy-light quark current qTQ with T being arbitrary Dirac matrices can be expanded in powers of l/mg as 



qVQ 



Here the contributions from both heavy quark and antiquark fields have been considered. 

According to above expansions for effective Lagrangian and effective current, one can write the matrix element in 
Eq. (2.1) as the following form in powers of l/mg, 



+ '-a^pF'^f')Qt\M,)+0{l/ml)}, (2.9) 



where Ab = ms — rat,, and F"'^ is the ghion field strength tensor. The effective heavy meson state \My) satisfies the 
heavy quark spin-flavor symmetry. Its normalization is 

{M,\Q+^^'Q+\M,) = 2kv^' (2.10) 

with the binding energy A = lim„jQ^oo Am being heavy flavor independent. 

It should be noted that the XjuiQ corrections in Eq.(2.9) include both contributions from the current expansion 
(2.8) and from the insertion of the effective Lagrangian (2.4). In Eqs.(2.8) and (2.9) the operator l/(w • D) arises 
from the contraction of effective heavy quark and antiquark fields [15,18]. In the v ■ A = Q gauge to be used in our 
calculation, this operator is tantamount to the heavy quark propagator. 

As can be seen in Eq.(2.9) that the l/mg order corrections to B — > plv transition are only attributed to one 
kinematic operator and one chromomagnetic operator. The heavy quark symmetry enables us to parameterize the 
matrix elements in HQEFT as 

(p(p,e*)|urQ+|M,) = -Tr\^{v,v)^M^l (2.11) 
(p(p,e*)|tZr-|+^7^iQ+|M,) = -Tr\^^{v,v)^M,l (2.12) 

{p(p,e*)\uT^^'-a^0F'^^Qt\M,) = -Tr[Qfiv,p)TP+'-a^0M.], (2.13) 
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where the pseudoscalar heavy meson spin wave function A4v = — is independent of the heavy quark flavor. 
^{v,p), Qi{v,p) and n"^{v,p) can be decomposed into Lorentz scalar functions, 

n{v,p)^Li{v p)/k* + L2{v ■ p){v e*) + [L^{v ■ p)fi* + Li{v p){v ■ e*)]p, 
fliiv,p) = 6Liiv ■ p)jk* + 5L2{v-p){v • e*) + [5L^{v ■ p)fi* + 5U{v p){v ■ t*)]p, 

^f{v,p) = {ri^ - fl") b*{Ri + R2P) + [v ■ e*)(i?3 + R4) 



+ (e*«7^ - e*^7") + Bsp + {e*^f - e*l^r) R7 + Rap 



+ ia' 



a/3 



/t{R9 + Riop) + [v ■ e*){Rii + R12P) 



with p'^ = p^^ /{v ■ p) and ^ = v-p = {m\ + — q^)/2mB- 
Eqs.(2.9)-(2.14) yield 

{p{p,e*)\uY{l-l'')b\B{pB)) = -2iy 



- {L', - L',){v ■ e*)r - L',{v ■ e*)v^ + 0(l/m^)} 



with 



R[ = -2Ri - 2R5 + R7- 3i?9 + (2i?2 - i?8)p^ 
iJ^ = ^2Rs - 2R5 - 3i?ii - (2i?4 + Rs)p 

/?3 = 2Ri — 2i?2 — 2Rq — i?7 + i?g — 3-Rio, 



s2 



R'a 



-2Rq — Rf — 2i?3 — 2i?4 — 3i?i2- 



(2.14) 



(2.15) 



Comparison between Eqs.(2.1) and (2.15) gives relations between the form factors Ai{i = 1, 2, 3), V and the universal 
wave functions, 



{L[{vp) + L'^{vp)} + - 



A2{q'^) = 2{mB +mp)^ 
Aaiq"^) =2{mB + mp)^ 



2m| 



2ms(u • p) 



ImBA^ L'^jv-p) L'^jv ■ p) - L'^jv ■ p) ^ ^ 



2m| 



V{q') 



I msA niB + nip , 



Ab 



2mB{v ■ p) 



(2.16) 



Ab mB{v-p) 

where the dots denote higher order l/mg contributions not to be taken into account in the following calculations. 



III. LIGHT CONE SUM RULES IN HQEFT 

For the derivation of the l/mg order corrections to B — > plv decay, we consider the two-point correlation function 

F'' = iJ d*xe-^(f--'"«'')-(p(p,e*)|r{u(0)7''(l - j')-^P+(dI + ^a„^F«/5)o+(0), 

Qtix)n'dix)}\0) (3.1) 

where Q^{x)ij°d{x) is the interpolating current for B meson. Inserting between the two currents in Eq.(3.1) a 
complete set of intermediate states with the B meson quantum number, one gets 
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itiqAb zAb — w L J 

+ / da ^^^' ^ + subtraction (3.2) 
J so 

with ^ = w • p and w = 2t> • fc, where k = pb — rnqv is the residual momentum of the bottom quark. The second term 
in (3.2) represents the higher resonance contributions. F is the decay constant of B meson at the leading order of 
l/mg expansion, defined by [17] 

{<d\qTQ+\B,) = ^Tr[TM,]. (3.3) 

In deep Euclidean region the correlator (3.1) can be calculated in effective field theory. The result can be written also 
as an integral over a theoretic spectral density, 

j ^g Pthii,s) ^ g^^^raction. (3.4) 
Jo S-UJ 

The standard treatment of sum rule is to assume the quark-hadron duality, and to equal the hadronic representation 
(3.2) and the theoretic one (3.4), which provides an equation: 

itiqAb 2Ab — CO ^ J 

+ r dsMlf) = r rf /*"^^'"^ + subtraction. (3.5) 

J so ^-<^ Jo S-UJ 

To ensure the reliability of sum rule estimates, one should enhance the importance of the groimd state contribution, 
suppress higher order nonperturbative contributions and remove the subtraction. These can be achieved by performing 
the Borel transformation 

4-). hm 

- -- n\ dhj 



4-)_!_ = e-Vr, B^e^^ = 6{X - (3.6) 



-ui/n = T 

to both sides of the equation (3.5). With using the formulae 

1 _„ /'r \,., .-. , , 1 

T' 
one gets 

2iF{5L'^€^''"''f^€lp,,Vfj - i{6L[ + 6L's)e*'' + i{5L'^ - 5L'^){v ■ e*)p^ + iSL'^iv ■ e*)?;^}e-2^B/T 

' dse-'/^p^^,s), (3.7) 

where the spectral density p(^, s) can be derived via double Borel transformations, 

p{^,s) = B[-y^^B!^^F^{^,u,). (3.8) 

In calculating the three-point function (3.1), one may represent the nonperturbative contributions embeded in the 
hadronic matrix element in terms of light cone wave functions. Among them are the two-particle distribution functions 
and the three-particle ones. However, if wo restrict our calculation to the lowest twist (twist 2) level, the three-particle 
functions do not contribute. As a result, the chromomagnetic operator in Eq.(3.1) could be neglected in the lowest 
twist approximation. The leading twist distribution functions are defined by [2,4-6] 

< p{p,e*)\u{0)a^M\0 > = -ifp{e;p. - <P^) C due'-'^-^Mu), 

Jo 

<p{p,e*)\u{0)-f^d{x)\0> =fpmpPi,^—^ [ due™f -^(^n («) 

P ' X Jo 

+ fprn,{e; - p/—^) f due'^^-^ g^l\u), 
<p{p,e*)\u{Q)^^^r,d{x)\Q> =\fpmpe^,^pe*"p''x^ j due'^^''^ g^l\u) (3.9) 
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with (t>A_,\\, and g^f'"'"' being functions with nonperturbative nature. 

Then the effective heavy quark fields Qy{xi)Qy{x2) can be contracted into a propagator of heavy quark, 
P+ dtS{xi — X2 — vt). In the lowest twist approximation, only the kinematic operator contributes to the l/mg 
order corrections to S — > plu decay. At v ■ A = gauge we used, the correlation function (3.1) simplifies as 



dl 



rfte-''=-(p(p,6*)|S(0)7^(l - l'')P+5{-y - vl) 



iv) 



+ '-a,,pF''^{y) + A„{y)A''{y)\p+5{y -x- vt)^'-d{xm 
d^x J d*y J dl J dte-'^-''5{-y-vl) - v"v \^^^^^^^ )6{y - x - vt) 



{p{p, e*)\u{0)^ [iV - 7" - ^"^(it^/35a^ + le^.a0V.)] d{x)\0) 



(3.10) 



with d"y-^ = d/{dya)- The final formula in (3.10) includes only the terms related to the two-particle distribution 
functions (3.9). 

Now the transition matrix element can be evaluated through the distribution functions defined in (3.9). Using Eqs. 
(3.6) and (3.9), the spectral function is found to be 



^/,m3(«^5f)" + />^i(«^0^)' 
1 



+ fpTUpiu^g'-l^y - /pmp(«20||)'] + {v ■ e*)v^ [2f^mlu4>^ - 2f^m%u'ct>,_y 



(3.11) 



where ' denotes derivative with respect to the variable it, while G^^\u) and $||(m) arc functions related to and (/>|| 

by ^G^^\u) = g''_^\u), ^$[|('u) = 4>\\{u). The detailed procedure in deriving Eq.(3.11)are similar to those in Refs. 
[10,11,13]. 



IV. NUMERICAL ANALYSIS 

and (^11 are the lowest twist distributions in the fraction of total momentum carried by the quark in transversely 

and longitudinally polarized mesons. They can be expanded in Gegenbauer polynomials Cn^'^{x) whose coefficients 
are renormalized multiplicatively. With the scale dependence explicitly, one has [4] 

</'±(||)(w,M) = 6u(l-«)[l+ Yl ai^\\\ti)C^J\2u-l)], 

n=2,4,--- 

ai^^H^^) = a;^(ll)(Mo)(^)(-""'*--o'"")/(^'^°\ (4.1) 

where Po = 11 — (2/3)n/, and 7JI, 7^ are the one loop anomalous dimensions [19,20]. The nonperturbative parameters 
and a| have been obtained in [4] with the values 

a^(lGeV) = 0.2±0.1, a|(lGeV) = 0.18 ± 0.10 (4.2) 
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and an'-^^^ = forn 7^ 2. 

The functions and g'f^ describe transverse polarizations of quarks in the longitudinaUy polarized mesons. As 
in Ref. [6] they are parameterized as 

9^l\u,n) = ^(1 + (2u - if) + ^af (m)(2u - 1)^ + (^4(m) + 56(/^))(3(2« - 1)^ - l) 
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15 



9±\u,n) = 6m(1 - u) 



^^^s{l^){3wX{f,) - uiifi))] (3 - 30(2u - If + 35(2« - 1)^) , 
- alip){2u - 1) + (ia|(M) + Iwil - ^Loiii,) + ^co^M)) 



(5{2u-lf - 1) 



(4.3) 



All the nonperturbative parameters in Eqs.(4.3) have been estimated in Ref. [6]. The asymptotic form of these 
factors and the renormalization scale dependence are given by perturbative QCD [21,22]. As in Refs. [11], the typical 
virtuality of the bottom quark 



ml - to2 « 2.4GeV, 



(4.4) 



is used for the energy scale for the current calculation. 

The values of the hadron quantities fp, fp, As, A and F have been extracted in the previous work (see, e.g., 
[4,17,23,24]). For consistency, here we use for them the same values as in Ref. [11], i.e.. 



fp± = (195 ± 7)MeV, fpo = (216 ± 5)MeV, = (160 ± 10)MeV, 
As w A = 0.53GcV, F = (0.30 ± 0.06)GeV^/2. 



(4.5) 



6Li as functions of T and sq can be derived from Eqs.(3.7) and (3.11). Fig.l shows the variation of 6Li as 
functions of the Borel parameter T at v ■ p = 2.5 GeV. The curves in each figure correspond to different values 
adopted for the threshold sq. 

The rule of LCSR method is to determine sq from the stability of relevant curves in the reliable region of T, where 
both the higher nonperturbative corrections and the contributions from excited and continuum states should not be 
large. In the current case, we focus on the region around T = 1.5 — 2GeV. As shown in Fig.l, SLi are found to be 
stable with respect to the Borel parameter T. In Ref. [11] the threshold sq = 2.1 ± 0.6GeV is adopted in evaluating 
the loading order wave functions Li. In calculating the decay width we will use for 6Li the same threshold values as 
those for the leading order wave function Li, i.e., sq ~ 2.1GeV. 

With (2.16), the form factors Ai, A2, A3 and V with including l/niQ order corrections can be calculated. It is 
convenient to represent each of these form factors in terms of three parameters as 



1 - aFq^/m% + bF{q^/m%) 



2 \2 ' 



(4.6) 



where F(g^) can be any one of Ai(g^), A2{q^), ^3(9^) and V{q'^). The parameters F{0), ap and bp presented in table 
1 are fitted from the the LCSR results at sq = 2.1GeV. Fig. 2 shows the form factors as functions of the momentum 
transfer, where the dashed curves are for the leading order results while the solid ones for the results with the l/mg 
order corrections included. 





m 


ap 


bp 




Ai 


0.26 


0.37 


-0.19 


LO 




0.27 


0.32 


-0.19 


NLO 


A2 


0.26 


1.11 


0.26 


LO 




0.26 


1.15 


0.30 


NLO 


A3 


-0.26 


1.12 


0.26 


LO 




-0.26 


1.11 


0.25 


NLO 


V 


0.32 


1.24 


0.29 


LO 




0.31 


1.26 


0.32 


NLO 



Table 1. Results of LCSR calculations up to leading (LO) and next leading order (NLO) in 
HQEFT. The leading order results are obtained in Ref. [10]. 
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The differential decay width of -B — > plu with the lepton mass neglected is 



X^/^^H^ + Hi + Hi) (4.7) 



and 



with the helicity amplitudes 

ms + rup 

Ho = -J-^{{ml -ml- q'){mB + mp)A,{q') ^A^iq')} (4.8) 

\={m\+ml^q^f -Ara^ml. (4.9) 

The total width oi B ^ plv can be obtained by integrating (4.7) over the whole accessible region of q^. We get 

r{B -> pliy) = (13.6 ± 4.0)|K6|'ps-\ (4.10) 

where the error results from the variation of the threshold energy Sq = 1.5 — 2.7GeV. 

The branching fraction of B° p'l+y is measured to be Br(B° p~l^v) = (2.6 ± 0.7) x lO"" [25]. This and the 
world average of the B° lifetime [25] Tbo = 1.536 ± 0.014 ps yields 

r(S° ^ p'l'^y) = (1.69 ± 0.47) X 10-*ps-^ (4.11) 

jV^bl is then extracted from Eqs.(4.10) and (4.11). It is 

IK&I = (3.53 ± 0.49 ± 0.52) x 10-^ (4.12) 

where the first and second errors correspond to the experimental and theoretical uncertainties, respectively. This 
value may be compared with the ones previously obtained [10,11,13,16,26]. Prom the exclusive semileptonic decays 
B 'K{p)ly, we then have 

1Kb I = (3.4 ± 0.5 ± 0.5) X 10"-'' {B nlv, LO) 

IK&I = (3.2 ± 0.5 ± 0.4) X 10"^ {B nlv, to NLO) 

IK&I = (3.7 ± 0.6 ± 0.7) X 10"^ {B -> plv, LO) 

\Vub\ = (3.5 ± 0.5 ± 0.5) X 10"^ {B plu, to NLO) 

\Vub\ = (3.5 ± 0.6 ± 0.1) X 10~^ {B inclusive semileptonic decays) 

As a summary, we have studied B plv decay up to the l/mq order corrections in HQEFT. In HQEFT, l/mg 
order corrections from the effective current and from effective Lagrangian are given by the same operator forms, 
which simplifies the structure of transition matrix elements. These l/mg order contributions have been calculated 
using light cone sum rules with considering the lowest twist distribution functions. Numerically, the 1/toq order 
wave functions give only corrections lower than 10% to the transition form factors. Similar to the B — » ttIv case, the 
correction indicates a slightly smaller value of the CKM matrix element \Vub\- The discussion concerning l/mg order 
corrections to B ^ plv decay in this paper is also applicable to other heavy to light vector meson decays. 
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Fig.l. Variation of l/mg order wave functions 5Li{i = 1,2,3,4) with respect to the 
Borel parameter T a,t = v ■ p = 2.5 GeV. The dashed, solid and dotted curves 
correspond to the thresholds so =1.5, 2.1 and 2.7 GeV respectively. 
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Fig.2. Form factors Ai(i = 1,2,3) and V obtained from light cone sum rules in 
HQEFT. The dashed curves are the leading order results in HQEFT [10,12], while 
the solid curves are the results with including 1 /mg order correction. 
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